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tH- ' We construct backreacted D3/D7 supergravity backgrounds which are dual 

to four-dimensional M = 1 and J\f = 2 supersymmetric Yang-Mills at large Nc 

with flavor quarks in the fundamental representation of SU(A^c)- We take into 

Qh! account the backreaction of D7-branes on either AdSs x S^ or AdSs x T^'^, or 

r~| , more generically on backgrounds where the space transverse to the D3-branes 

is Kahler. The construction of the backreacted geometry splits into two stages. 
• 1—1 . 
^ ■ First we determine the modification of the six-dimensional space transverse to 

C^ ■ the D3 due to the D7, and then we compute the warp factor due to the D3. 

The N = 2 background corresponds to placing a single stack of Nf D7-branes 

in AdSs x S^ . Here the Kahler potential is known exactly, while the warp factor 

is obtained in certain limits as a perturbative expansion. By placing another 

D7' probe in the backreacted D3/D7 background, we derive the effect of the D7- 

branes on the spectrum of the scalar fluctuations to first order in Nf. The two 

systems with N = 1 supersymmetry that we discuss are D3/D7/D7' and D3/D7 

on the conifold. In both cases, the Kahler potential is obtained perturbatively 

in the number of D7-branes. We provide all the ingredients necessary for the 

computation of each term in the expansion, and in each case give the first few 

terms explicitly. Finally, we comment on some aspects of the dual gauge theories. 



1 Introduction 

In the framework of the gauge/gravity correspondence [1] a complete holographic de- 
scription of QCD is still beyond our current reach. Nevertheless, substantial progress 
has been achieved in the holographic understanding of pure A/" = 1 super Yang-Mills 
(SYM) at large Nc [2,3]. The holographic understanding of flavor dynamics at the 
quantitative level has developed slower, in great measure due to the lack of a fully back- 
reacted supergravity background describing gauge theories with fundamental flavor. In 
this paper we address the construction of fully backreacted supergravity backgrounds 
dual to A/" = 1 SYM with fundamental flavor. 

Recently, there has been much attention devoted to flavor dynamics in the probe 
approximation limit [4-16]. However, a signiflcant limitation of the probe approxima- 
tion is that in that regime the number of flavors is necessarily much smaller than the 
number of colors: Nf/N^. ^ 1. Still, even in the probe approximation, there have been 
some quite interesting results [17,18]. 

Why go beyond the probe approximation? One of the major stepping stones in a 
holographic understanding of QCD is a control of flavor dynamics. For example, dy- 
namical mesons and baryons play a central role in obtaining a complete understanding 
of the QCD phase diagram. In particular, in the large Nc limit the critical temper- 
ature should behave very differently in the presence of baryons. On the other hand, 
it is known that baryons play a subleading role at large N^., which leads to the ques- 
tion of how precisely they change the critical temperature. On a more formal level, 
the quantum moduli space of SU(A'"c) supersymmetric gauge theories with Nf flavors 
crucially depends on the relative value of Nc and A'^^. There is a host of phenomena in 
flavor dynamics of supersymmetric theories, such as the conformal window, that can 
be explored only in the limit of comparable Nc and Nf. Having Nf <^ Nc excludes, 
in principle, the possibility of a holographic description of such phenomena. Clearly, 
unless a fully backreacted solution becomes available, a quantitative understanding 
of several physical questions of gauge theories will remain elusive in the holographic 
approach. 

In this paper we address the question of backreaction in the context of supergravity 
backgrounds dual to A/" = 1 supersymmetric gauge theories. While the general class of 
supergravity backgrounds dual to A/" = 1 gauge theories is wide, we concentrate only 
on backgrounds containing D3 and D7-branes. In this case the backreacted metric of 



the space transverse to the stack of D3-branes remains Kahler. For Kahler metrics a 
remarkable simphfication of the Einstein equations takes place. In fact, solving the 
Einstein equations is equivalent to solving an equation that in many cases becomes 
a standard Monge-Ampere equation. Therefore, the problem of finding backreacted 
solutions splits into two main steps. In the first one, after finding a suitable Ansatz 
for the dilaton, we determine the Kahler potential of the resulting six-dimensional 
background transverse to the D3 worldvolume. In the second step we find the warp 
factor due to the presence of the D3-branes. 

The paper is organized as follows. In section 2 we review the construction of the 
D3/D7 system to introduce notation and comment on some of the choices that we make 
further in the paper. Section 3 contains a description of the D3/D7/D7' system, that 
is, two perpendicular stacks of the D7-branes that are parallel to the D3. This system 
has A/" = 1 supersymmetry, and we discuss it in some detail. Section 4 presents a probe 
calculation of the scalar (meson) spectrum of the A/" = 1 gauge theory. We provide the 
spectrum as an expansion in the number of D7-branes of the D3/D7 geometry. We then 
turn to the conifold picture by first reviewing the structure of the conifold as a Kahler 
manifold in section 5, and then describing a solution with D7-branes holomorphically 
embedded in the conifold in section 6. Section 7 contains some comments on the gauge 
theory side. We conclude in section 8 by enumerating a number of open problems. 
Appendix A contains some useful results pertaining to properties of the Monge-Ampere 
equation. 

2 D3/D7 Systems 

In this section we review the basic setup for examining D3/D7 systems. The D3- 
brane geometry by itself breaks the 10-dimensional Poincare invariance into manifest 
SO (3, 1) X SO (6). This may be represented by a metric of the form 

dsjo = h-^'\ym)dxl + h'/\ym)gmndy"'dy'', (2.1) 

where x^ (yU = 0, . . . , 3) are fiat coordinates on the longitudinal spacetime, and y™" {m = 
4, . . . , 9) are coordinates transverse to the D3-brane. The D7-branes are embedded in 
this 10 dimensional space such that they share the four longitudinal directions with the 
D3-branes while wrapping or spanning four out of the six transverse directions. This 
embedding may be represented pictorially as 
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for the case of a single stack of D7-branes. 

D7-brane configurations have been effectively studied in [19] in the context of cosmic 
strings, and D3/D7 systems have been examined in [20,21]. The starting point for 
exploration of D3/D7 geometries preserving some fraction of supersymmetries is the set 
of Killing spinor equations obtained from the JIB supersymmetry transformations [22] 
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are the SL(2,M) scalar kinetic and composite connection terms for the dilaton/axion 
pair T = Ti + iT2. 

The Killing spinor equations simplify in the absence of G^ flux to yield 
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Corresponding to the 4 + 6 coordinate split of (2.1), we decompose the Dirac matrices 
according to 7^ = (F^ 1), 7™ = (F^ (g) F™"), where F represents either the S0(1, 3) or 
SO (6) Dirac matrices as appropriate. This allows us to deflne the chirality matrices 
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The complex JIB spinor e has deflnite chirality, 76 

„5 



e. Hence the four and six- 
dimensional chiralities are correlated, so that l^u)^ = iJq)^ = se where s = ±1. In 
addition, the coordinate split also allows us to write the most general flve form consis- 
tent with four-dimensional Poincare symmetry and self-duality: 
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where F^ is a function only of the y™ coordinates, and is zero for M = . . . 3. To insure 
Poincare symmetry we must also take r to be a function only of the y™' coordinates. 
Using the above, the Killing equations, (2.4), can now be written as 

P^(l®r™)e* = 0, 
9^e-('^9„log(/i) + iF„Vr^®ne = 0, 

V^^e - '^F„,e + r-dn log(/i) + '-FA (1 ® r„")e - '-Q^e = 0, (2.7) 

where V*-^-* is covariant with respect to the metric on the six dimensional space trans- 
verse to the D3 branes. Taking e to be only a function of y^, the second equation can 
immediately be solved by taking 'yf^\^ = e (so that s = 1) and at the same time setting 

F^ = -]dm\og{h) = -^. (2.8) 

4 An 

This is a remarkable result, as it demonstrates that the usual relation between the D3- 
brane warp factor h{y"^) and the five-form remains unchanged, even in the presence of 
D7-branes. Note that the s = 1 condition on Killing spinors yields the familiar result 
that the D3-branes by themselves are half-BPS objects. 

Taking the warp factor as above, the Killing spinor equations (2.7) reduce to 

7"^P^ry* = 0, (V„ - '-QJv = 0, (2.9) 

where t] = h^/^e. These equations may be used to determine the structure of the 
D7-branes. Since the second equation above implies that r] is covariantly constant, the 
connection Vm is in U(3), and the transverse space to the D3-branes is complex Kahler. 
The first of the above equations is a familiar one for D7-branes, namely {'j'^dm.T)'!]* = 0, 
and can be satisfied by taking r to be holomorphic on the complex 3-fold. This yields 
the further conditions 'y"^ri* = on the Killing spinors whenever S^r ^ (for the 
generic case), where z^ {in = 1,2,3) is now a complex coordinate. As a result, the 
D3/D7 system preserves 1/4, 1/8 and 1/8 of the JIB supersymmetries for r depending 
holomorphically on one, two and three of the complex Zm coordinates, respectively. 
Note that there is no further reduction of supersymmetry for the last case, as the 
product of the three independent D7 projections gives precisely the D3 projection. 

The second equation in (2.9) has a useful integrability condition, Rmli = Pm^n + 
P*nPn- This has a particularly nice form for Kahler geometry, namely Rmn = PmPn^ 



where m and n denote complex indices. Taking T{zm) to be holomorphic and using Pm 
from (2.3) results in 

-9„9^1og(det((7(^))) = ~d^dn\og{T - r*), (2.10) 

which may be integrated to give 

det{dmdnJ^) = nnim{T). (2.11) 

Here JF is the Kahler potential and Q is an arbitrary holomorphic function accounting 
for the integration factors. The same Q is used for the holomorphic and antiholomorphic 
integration functions to insure that the Kahler potential is real. In cases with 'radial' 
symmetry, Eq. (2.11) can be reduced to a real Monge Ampere equation (see Appendix A 
for further comments on this equation). 

Finally, returning to the D3-brane warp factor, the integrability condition combined 
with the Einstein equations gives the result that the function h is harmonic 

n(6)/i = -i27T)^NAetig^'^y'^'6' {y^ - y^,) . (2.12) 

As a result, solving for the D3/D7 geometry can be reduced to several steps. First, 
we take a (noncompact) Calabi-Yau background geometry, and then turn on the D7- 
branes by choosing a holomorphic r with appropriate monodromies. Then we must 
solve the complex Monge Ampere equation (2.11) for the Kahler potential, yielding 
the geometry of the six-dimensional transverse space. Finally, to include the D3- 
branes, we must solve the harmonic function equation (2.12). This gives the warp 
factor and corresponding five-form flux through (2.8). Unlike for most intersecting 
brane configurations, at least in principle this procedure yields fully localized D3/D7 
solutions, although the resulting equations are often quite intractable in practice. Note 
that the possibility of constructing localized intersecting solutions for the D3/D7 system 
is a feature shared with the D2/D6 system [23]. 

2.1 D3 and a single stack of D7 in fiat space 

We now consider the simplest case where the six-dimensional transverse space to the 
D3-branes is flat before turning on any D7-branes. In particular 

C^Se, background = dZidZi + dZ2dZ2 + dz^dz^, .^flat, background = ^1^1 +^2^2 + ^3^3- (2.13) 



In the presence of D7-branes, the Kahler function may be obtained by solving (2.11). 
However, before doing so, we must understand the form of the dilaton/axion field, 

t{zi,Z2,Z3). 

For a single stack of D7-branes transverse to Z3, we may follow the procedure given 
in [19], which we outline here. If one requires modular invariance for r, holomorphicity 
and the proper monodromy around the D7 branes, the basic solution can be written 
in terms of the modular invariant j-function as [19,21] 

J{r) = fizs), (2.14) 

where for now / is an arbitrary holomorphic function. In [19] it was observed that poles 
in / correspond to D7-branes, in the sense of decompactification at the cores of the 
branes (as will be indicated below). Solving (2.11) with this form of the dilaton/axion 
is quite involved; hence we only consider the neighborhood close to one of the poles 

The role of the j'-function in the above expression is to provide a map from the funda- 
mental domain into the whole complex plane which is then related to Nf coverings of 
the zs-plane. 

We have limited ourselves to z^ close to zq which requires that r2 is large. In this 
case the j function can be approximated as exp(— 27rir), so that the above equation 
becomes 

r(^3) = ^ (^ - Sl°g(^3)) , (^3 « 1), (2.16) 

where we have shifted z^ to remove Zq. This is known as the "decompactification limit" 
from F-theory because it corresponds to the modulus of the compactified torus going 
to infinity. We will work in this limit for the remainder of the paper. (Note that 
Eq. (2.16) is no longer appropriate for Z3 > 1, where instead r approaches a constant 
value.) It is also natural to consider scalings of z^ so that b may be removed as well. 
Two scalings that are natural are b = 1 and b = e^'^, making the constant log(6)/27r 
either or 1, respectively. 

In [19], Q is fixed by considering modular invariance of the metric, imposing regu- 
larity of the metric, and noting that (733 is proportional to T2. This gives 



933 



T2\v{T)\'\zs\-''f/'. (2.17) 



Near the core of the D7-brane, where T2 is large, this has the simple form 

933^^2- (2.18) 

In particular, the additional factor |^(T)|'^|2;3|~^-f'^ in the metric simply approaches a 
constant near the cores of the D7-branes. 

We now turn to the question of the Kahler potential. We have chosen to stack the 
D7-branes together so that they are extended along zi and Z2, and are transverse in 
Z3. Making use of the background symmetry, and noting that T2 is a function of the 
real combination Z3Z3 suggests that we take 

J^D7 = Z1Z1 + Z2Z2 + fiz^h)- (2.19) 

The background Kahler potential is obtained from (2.11) by taking Imr = 1 and Q = 1, 
so that when / = z^z^ we simply get flat space. We now take b = e^'^ and turn on 
non-zero Nf, so that (2.11) becomes 

9393/(^3^3) = iz,z,)f' + f=(l-^ \ogiz,z,)] . (2.20) 

By making a change of variables, ys = log{z3Z3), we arrive at 

/"(l/s) = (1 - ^1/3) ey-\ (2.21) 

This may be straightforwardly integrated to obtain the Kahler function for a single 
stack of D7-branes located at Z3 = 

J^D7 = ZiZi + Z2Z2 + Z3Z3 - -^{z3Z3){l0g{ZsZ3) - 2) + Ci logizsZ^) + C2. (2.22) 

An 

The constants Ci and C2 do not affect the metric; we will therefore set them to zero (we 
may always add a holomorphic plus antiholomorphic function to a Kahler potential). 
Of course, the transverse space metric is given by 

dsl = dzidzi + dz2dz2 + ( 1 — 7^ log -23^3 ) dz^dz^, (2.23) 

and this could have been obtained directly from (2.18) without even having to solve for 
the Kahler function. Nevertheless, this exercise in obtaining JF by solving (2.11) will 
prove essential when considering multiple stacks of D7-branes as well as D7-branes in 
other backgrounds. 
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2.2 The warp factor 

Now that we have found the D7-brane contribution to the transverse metric, (2.23), we 
may turn to the D3-brane warp factor. To determine the warp factor, we must solve 
the transverse Laplacian (2.12). For a single stack of D7-branes, we first rewrite the 
metric (2.23) in the generic form 

dsl = dzidzi + dz2dz2 + e^^^^'^^'Mz^dz^, (2.24) 

in which case the transverse Laplacian is more explicitly 

n(6)/i = {didi + 828-2 + e-''8^8i)h. (2.25) 

This expression for the warp factor was examined in [24] , where it was solved to a first 
order approximation away from the D7-branes, where the metric is approximately fiat, 
with a deficit angle. 

In the absence of D7-branes, it is natural to place the stack of D3-branes at the 
origin, Zm = 0, in which case the S0(6) R-symmetry is preserved, so that h = L^/R^ 
where R is the radial coordinate, R^ = |-2iP + |-22p + ksP- In orienting the D7-branes 
to lie fully in the zi and Z2-planes, we preserve a natural SO (4) x SO (2) C SO (6) 
subgroup. Note, however, that separating the D3 and D7-branes in the zs-plane will 
break the SO (2). With this symmetry in mind, we introduce real coordinates 

r' = kir + k2r, ^3=P3e^'^^ (2.26) 

so that 

°{6)/^ = (^dy8r + e-*(^^'^«)(-5p3P35,3 + \8lj) Mr,P3,03). (2.27) 

V^ P3 pi J 

If the D7-branes are separated by a distance d along the real zs-axis, the metric function 
e* has the form 



e^ = ( 1 - -^ log(p^ - 2psd cos 03 + d')\. (2.28) 

Although no straightforward solutions to (2.27) appear to exist in the general case, 
we may find approximate solutions in appropriate limits. For example, when d = 
the D3-branes overlap the D7-branes, and the remaining SO (2) symmetry is restored. 
In this case, the Laplacian reduces to one with only two radial variables, p^ and r. 
Furthermore, we may consider solving for the warp factor near the core of the D7- 
branes, where e* ^ 1. This is the decompactification region of the complex zs-plane, 
and corresponds to the region where the backreaction is strong. 



It can be checked that the equation^ 2.27 allows separation of variables in the form: 

h = f{r)g{p,<j>)Y\S'), (2.29) 

where Y\S^) are spherical harmonics on S^ 

V*Vi F' = -/(/ + 2) F'. (2.30) 

The equation for the function /(r) can be solved in terms of Bessel functions. The 
equation for g can be further separated in variables, after Fourier expanding in an 
appropriate angular variable, resulting in an equation for the radial variable along. 
The latter is a more involved second order differential equation which we will not 
discuss here in detail. 



3 D3 and two stacks of D7-branes 

The D3 geometry in the presence of a single stack of D7-branes is perhaps the most 
straightforward to investigate, as there is little ambiguity in how the D7-branes are 
turned on. However, we are generally more interested in configurations admitting 
A/" = 1 duals with flavor. This may be achieved by either starting with a transverse 
space to the D3-branes with reduced supersymmetry or by turning on multiple stacks 
of D7-branes. In this section we focus on the latter possibility by turning on two stacks 
of mutually perpendicular D7-branes. 

In particular, we can construct a solution representing two stacks of D7-branes by 
placing, as before, N3 D7-branes perpendicular to the zs-plane and another set of N2 
D7-branes transverse to the 2;2-plane. This configuration may be represented as follows 
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Note that turning on a third stack of D7-branes transverse to Zi would not further 
reduce any supersymmetry. However for simplicity we avoid this case as well as other 
possible non-perpendicular configurations. 



^We are thankful to A. Hashimoto and M. Mahato for discussions relevant to this topic. 



In order to have two stacks of D7-branes, located at Z2 ^ and 2:3 ^ 0, we take 
t{zi, Z2, Z3) to have the appropriate monodromies around the respective cycles. More 
precisely, near the D7-branes we have 

r = .(l-§log.2-§log.3). (3.1) 

As a result, it is natural to take an Ansatz for the Kahler function in the form 

^D7/DT = ZiZi + f{z-2Z-2, Z^Z^). (3.2) 

Since both J?^d7/D7' and Imr depend only on the magnitudes of the complex coordinates, 
we make a similar change of variables as in the single stack case 

1/2 = log(2;2^2), 1/3 = log(^3^3), (3.3) 

so that the equation determining the Kahler potential, (2.11), becomes 

/22/33-/23^=(l-fl/.-§i/3)e---, (3.4) 

where subscripts on / denote differentiation with respect to the appropriate 1/2 or ?/3 
variable. This is a non-linear second order differential equation, and to our knowledge 
there are no general methods for obtaining a solution. However, as we explain in 
Appendix A, Eq. (3.4) is a Monge-Ampere equation in the two real variables 1/2 and 
1/3. For such equations there are general existence theorems for solutions, given some 
fairly general boundary conditions (see Appendix A). 

A crucial role in understanding the solutions of (3.4) is played by the choice of 
boundary conditions. In the following subsection we present a solution that is pertur- 
bative in the number of D7-branes. The boundary conditions that we impose imply 
that to lowest order the solution behave like two independent superposed stacks of 
D7-branes. 

3.1 Perturbative solutions for two stacks of D7-branes 

In order to develop a series solution to (3.4) where both A^2 and A^3 are small, we 
introduce an expansion parameter e and define 

N2 N, 

n2 = e-, n, = e-. (3.5) 
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We have also absorbed the factors of An to simphfy the subsequent expressions. Note 
that this expansion presupposes that both N2 and A^3 are of the same order. In this 
case, (3.4) takes the form 

/22/33 - /23' = ey^+y^il - e(n2|/2 + nat/s)). (3.6) 

We now proceed to look for a perturbative solution of the form 

00 
/(l/2,l/3)=/(°) + E^"/^"^- (3.7) 

n=l 

The zeroth order solution corresponds to flat space with no D7-branes, and the 
corresponding part of the Kahler potential is simply given by 

f(o) = ey^ + ey\ (3.8) 

At the first order, we find 

e-''fS + e-^'fS = -^21/2 - nm- (3.9) 

In general, solutions to this equation are only determined up to choice of boundary 
conditions. However, at the linear order, this can be viewed as a sum of two equations — 
one for each D7 stack. Hence we choose a simple solution of the form 

fW = -n^ey^y^ - 2) - nse^^d/g - 2). (3.10) 

This may be viewed as defining our boundary conditions (at first order) corresponding 
to the desired two stack configuration. 

We now demonstrate how this series expansion may be developed to higher orders. 
Although we have started with a non-linear equation, the perturbative expansion yields 
linear equations at each order. In particular, for n > 2, /'^"Ms determined by solving 

ra-l 

e-y^f!£^ + e-^3;W = _e-..-.3 ^ (/W/(3-) - /»4-)) , (n > 2). (3.11) 

Note that the linear operator L = e~^^(9| + e~^^(9| separates, allowing a solution to 
be obtained through separation of variables. However we will not pursue this method 
here, but will instead develop the next few orders by inspection. The second order 
equation becomes 

e-'^A? + e-^'fif = -(^2^3)1/2 1/3, (3.12) 

11 



demonstrating that 1/2 and t/s become mixed at this order. This equation has a solution 
of the form 

/(2) = ae^^l/3(l/2-2) + 6e^3l/2(t/3-2), a + b = -n2n3. (3.13) 

Of course, this solution is not unique. However, it maintains the structure of the lower 
order terms, namely / ~ 6^^/2(1/2,1/3) + 6^^/3(1/2! l/a)- In order to maintain symmetry 
between the two stacks, we take a = b = — {71211^) /2. Similarly, we find a solution to 
third order 

which maintains interchange symmetry between 1/2 ^^ Vz (with corresponding n2 ^^ n^)- 
Thus the above mentioned structure of / as a sum of two independent terms persists 
to third order. It is worth remarking, however, that the solution necessarily has a 
different structure at fourth order in perturbation theory. The reason is that precisely 
at this order the term /23^ starts contributing, and therefore breaks the symmetry of 
the first three orders. 

Collecting the above terms in the expansion, we find that the Kahler potential takes 
the form 



^D7/D7' = Z1Z1+ Z2Z2 + ^3^3 
N2 

N2N3 



N2 _ _ N3 _ 

--—Z2Z2 (logZ2^2 - 2) - -;— ^3^3 (log 23^3 - 2) 
An Air 



327r2 


1287r3 



Z2Z2{\ogZ2Z2 - 2) log Z3Z3 + 2:3^3 (log ^3^3 - 2) logZ2^2 
^2^2 log 2:3^3 (log^ ^2^2 - 4 \0gZ2Z2 + 6) 

1287r3 ^^^^ ^°^ ^^^^ ^'^'^ ^3^3-4 log 2:3^3 + 6)h . (3.15) 

Given the Kahler potential, we can easily read off the metric. We present a few entries; 
the rest can be recovered using the symmetries under Z2 and Z3: 



iV2, _ iV2iV3, _ , _ iV|iV3, 2 
log Z2Z2 7- log ^2-22 log Z-iZ-i '"" 



fl'22 = 1 - :;— log ^2^2 - ^^ ^ log ^2^2 log 23^3 - -,1^ o log^ ^2^2 log ^3^3 + 



932 = 577^ — (log ^2 ^2-1) + — (log 2:3^3-1^ 

327r^ \zs Z2 

2 3 A (iog2 ^2^2 - 2 log ^2^2 + 2) - -f-4 ^ (log2 ^3^3 - 2 log ^3^3 + 2) 



1287r3 Z3 ^ ^ 1287r3 Z2 

+ ■■■. (3.16) 
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We see that, to first order in the number of D7-branes, the metric looks hke a simple 
superposition of two stacks of D7. This validates our implicit choice of boundary 
conditions imposed above. However, already at second order we see a mixing between 
the Z2 and Z3 planes. Note that (732 signals a highly nontrivial dependence on the 
coordinates Z2 and z^ through the presence of ratios of complex coordinates. 

3.2 A perturbation in the relative number of D7-branes 

While in the above we have focused on N2 and A^3 both of the same order, we may 
alternatively consider the case where one of the stacks of D7-branes is much smaller 
than the other. To be concrete, we take N2 '^ N3, and consider an expansion in terms 
of the small parameter N2/N^. In this case, we define 

N2 Ns .__. 

^^ = '^^ -3 = ^, (3.17) 

so that the determining equation, (3.4), becomes 

/22/33 - /23' = e^^+^^((l - n^ys) - 6^2^/2). (3.18) 

Clearly this equation may be expanded in the same manner as (3.6). However, the 
source terms will be rearranged at the first two orders. 

At zeroth order, we take the exact solution for a single stack of A3 coincident 
D3-branes 

fio) = ey^+ey'{l-n,{y,-2)). (3.19) 

At first order, we find that f'^^^ must satisfy 

e-y'^{l - nm)^^ + e-y^f^l^ = -7221/2. (3.20) 

This has a particularly simple solution of the form 

/W = -n2ey'y2. (3.21) 

For higher orders, the perturbative expansion has the form 

{^-"^'^ + r^^) /'"> - -S^ g (/^f^-" - ^SVr") , (" > 2). 

(3.22) 
which is again amenable to separation of variables. 
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Collecting the terms in the Kahler potential computed above, we obtain 

N3 _ _ N2 _ _ 

T = ZiZi + 2:22:2 + 2:32:3 - -— Z3Z3(logZ3Z3 - 2) - -— Z3Z3logZ2^2 + OiN^). 

4:71 4:71 

This Kahler potential results in the following metric: 

ds"^ = dzidzi + dz2dz2 + ( 1 log 2:3^3 \0gz2Z2 ] dz^dz^ 

\ 47r 47r / 

-^ (-dz3dz2 + -dz2dz3] + OiNi). (3.23) 

47r \Z2 Z2 J 

Note that the main effect of turning on a second stack of D7-branes transverse to 2:2 
is to change (733 by an expected amount proportional to log 2:2^2- What is slightly 
unexpected, however, is the appearance of non-diagonal terms (723 at first order in the 
N2 perturbation. Nevertheless, this is consistent with the expansion of the previous 
subsection. 



4 Probe analysis 

In the previous section, we have examined the Kahler structure of the background 
with two perpendicular stacks of D7-branes. To proceed, we would now in principle 
need to solve for the D3-brane warp factor, much as we did for the single stack case 
in section 2.2. However, we may in fact already extract some A/" = 1 fiavor physics 
without full knowledge of the warp factor by considering a perpendicular D7' probe in 
the backreacted D3/D7 background. 

Let us first reemphasize, however, that taking the F-theory expression for the dila- 
ton, (2.15), implies that the dilaton/axion has a very different behavior in two regions 
of interest: near and far from the D7-branes. Near a stack of D7-branes we have an 
expression of the form (2.16). However, far away from the stack (that is, for large 
values of z) the solution of (2.15) is roughly of the form 

^ = ^0 + ^. (4.1) 

The above expression can be obtained from (2.15) by considering an expansion near 
\z\ -^ 00. Note that the radius of convergence of the above expression excludes z = 0. 
The leading term is nothing but a constant dilaton thus allowing us to recover fiat 
space in this limit. 



14 



The general solution to (2.11) in the case of one stack of D7-branes implies that the 
metric entry Qzz is simply proportional to T2. The arbitrary function Vt of (2.11) was 
fixed in [19] by imposing modular invariance and regularity of the metric, leading to 



Qz 



r,\r^{r)t\z\-''^"- 



(4.2) 



In the large \z\ limit that we are interested in in this section, the metric becomes that 
of flat space with a deficit angle. Namely, for large \z\ the dilaton is constant to first 
order and therefore we have for the relevant part of the metric (defining z = re**^) 



ds 



r-'^f/'idr' + rW) 



1 



[i-^y 



N 



dp' + Al--^ 



f\2. 



(4.3) 



where p = r^-^s/^"^. Thus the effect of the D7 is to create a deficit angle at asymptotic 
infinity. Due to this deficit angle we are forced to consider Nj < 12 (where p = logr 
must be taken for Nf = 12)^. 

To be specific, we take the orientation of the probe D7' brane as indicated in the 
table below 








1 


2 


3 


4 


5 


6 


'~7 

1 


8 


9 


D7 


— 


— 


— 


— 


— 


— 


— 


— 






D7' 


— 


— 


— 


— 






— 


— 


— 


— 



The D3/D7 background is 

ds'^ = h~^''^dsl + h^''^{dzidzi + dz2dz2 + e^dz^dzs] 



a 



0123 



-h~\ 



n(6)/i = -(27r)^iV,e-^5^(^). 



(4.4) 



In the static gauge, with an Ansatz such that angular coordinate 0i in the 2;i-plane is 
constant and the radial coordinate has a profile pi = pi{pz)i the pull-back metric on 
the D7' probe becomes 



*gabda''da'' = h~^'^dsl + h^'^{dz2dz2 + (e* + pi^)dpl + e^p^^^s^), 
where pi = -j^Pi- The action of the D7' probe 



(4.5) 



S=-pj j d^ae-^^-det{*g + 2T,a'¥)+pj j *C^,^ + ^S^ j *C^^^¥^¥ (4.6) 

■^Actually, two asymptotic geometries with Nf ~ 12 may be joined to provide the F-theory fibration 
of K3 over a compact base. However, we do not consider this case. 
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becomes (up to some infinite volume factor 2'n-Vol{R^)) 

S=-firl rfp3e-*e*/V3A/e* + Pi'- (4-7) 

Tlie equation of motion for pi determines tlie profile of the probe brane 

^ I e-^e-^/^ ^i^ I = 0, (4.8) 

where in the asymptotic region e* = P3 ^ and e* = const^. The most general 
solution in the asymptotic region is given by 



Pi = I dp3 ' 



Ps ~ "-"iPs 



/CT 



N-12 



({N - 12) logp3 - 12 log (1 + ^J^~C^i^^l^\\ + C, (4.9) 



where Ci,C2 are integration constants. 

Since pi = const, is a solution, the probe can be placed at a fixed position in the 
Zi plane and therefore the probe has a flat profile. A consequence of this fact will be 
the absence of a quark condensate, i.e. (QQ) = 0. 

4.1 The spectrum of the scalar fields 

By expanding the probe action in fiuctuations around the fiat profile solution pi = 
d, 01 = const., which is a solution even in the presence of a running dilaton, and 
identifying the transverse fiuctuations X, Y with the scalar fields of the dual gauge 
theory, we extract their action (to leading order) as 

S=-fij f v/dit^e-*(l + *g''\d,XdtX + d^YdbY)), (4.10) 



where 



Vdet *g = p2p3e'^, e* = 1 - -^log(p3), 

Zn 

e-* = l-^log(p3). (4.11) 



^From T — Too + -, — Tvre *^/'^3 we extract the dilaton e * = Inir = IniToo §- siniNfch^). 

(p'i) \ ^ ^ ^ P. . 

Furthermore, upon integration over ^3, we remain only with the constant dilaton factor in the effective 

action of the D7' probe. 
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The equation of motion for either of the two scalar fields (collectively denoted by X) 
is 

"!..„,,. 1 



hri'^'d^d^X + 



—dp^{p2dp^X) + ^d.X 



,P2 



+e" 



P2 

-dp,ipsdp,x) + -^dl.x - dp,^dp,x 

.P3 Pz 



0, (4.12) 



where we recall that the warp factor h obeys 



{d,di + d2d-2)h + e-* (-dp,p:,dp, + \dl\ h = -{2n)^N,e-''6%z). (4.13) 

\P3 Ps / 

When substituting the warp factor into (4.12) we evaluate it at a fixed position in the 
2;i-plane, namely h = h{pi = d,(f)i = const., p2, Ps). Even though the warp function 
equation cannot be solved exactly, we can attempt to solve it perturbatively in Nf, 
the number of D7-branes. With the change of variable e^/^ps = r^, the metric in the 
Z3-plane becomes 



P3 



rs ( 1 + — logrg 



dr^ + ridt 






and the corresponding Laplacian can be approximated by 

iV/\ 1 _ _ 1 



(4.14) 



e ^8383 



1 - -- —drsirsdr,) + -^d, 



27r / rg 



2 (ps' 



(4.15) 



The effect of the D7-branes on the background, to first order, is to generate a conical 
deficit in the otherwise fiat z^ directions. 

One way to proceed is to take full advantage of the observation that, to first order 
in Nf, the geometry of the six-dimensional space transverse to the D3-branes is (con- 
formally) fiat, with a conical deficit in the zs-plane. The warp factor is then simply 
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the Green's function [25] G{z, 0) associated with this six-dimensional geometry: 

dsl ^ dzidzi + dz2dz2 + dr^ + rldip'^, 

in ~ (Y) -|- 1 — ^ 27171 = v? + 2z/7r, 
V 47r / 

G{rlpl + pl^) -- 



pI + pi + '^i . 



sin((/9 — z/tt) 2(1 + cos((/9 — z/tt)) 

=arctanh- 



^2cos(2v9-2z/7r) -2 y/2 cos(2v9 - 2z/7r) -2 

sin ( (/9 + z/tt) 2(1 + cos((/9 + z/tt)) 

=arctanh- 



2 cos(2^ + 2z/7r) - 2 ^2 cos(2(^ + 2z/7r) - 2 



'4r^/ 2 „2 , .2n_ r4/i ^/^ 1 



. « L^G(.3^p^ + </^)«L*^l-^j^^^^;j^. (416) 

Thus to first order in Nf the effect of the D7-branes on the warp factor is a mere 
rescahng with (1 — Nf/2). This will carry through when investigating the meson 
spectrum. 

As we discussed earlier, the equation of motion for the two scalar fields correspond- 
ing to fiuctuations in the transverse directions Zi, 

Xir:,,ip,pl + d^)r^Xir:,,ip + 2u7r,pl + d^), (4.17) 

can be written as 

{n^ + M^h)X = 0. (4.18) 

where we have Fourier expanded the scalar fiuctuations in the directions common with 
the D3-branes, X = J d'^ke^^''^'^X{k), and we have identified the Casimir invariant k"^ 
with the four- dimensional mass k"^ = — M^. We continue to expand X into Fourier 
modes, according to the identifications of the wedge geometry 

X = 5^X„,„(p2,r3)e*'"^^e*'^"'^. (4.19) 

m,n 

By substituting into the equation of motion, we find 



p2 r^ Mir^ P2 Tg 



—dp,p2dp, + -dr.r^dr, + ^dr, -'-^-^-M^h] X™,„ = 0. (4.20) 



The next step is to make yet another change of variables, r^ = i^cosa, p2 = -Rsina, 



and to further expand Xm,n{R, «) into the eigenfunctions of the differential operator 

1 Ti V TTL 

C? = daisin a COS ada) 7. 5 — , 

sm a cos a cos^ a sin a 

C^ia) = Aa(a), 

X{R,a) = X{R)Uc^), (4.21) 

where A = — /(/ + 2) and / = \m\ + |n|z/. We have thus reduced the equation of motion 
of the scalar fluctuations to an ordinary differential equation for X{R): 

■^dnR'dR + ^^ + ^dn + M^h\x{R) = 0. (4.22) 

The square normalizable solution is 

AT, 1 

X{R) = R-^-^+\R^ + d^)-\Fi{-a,- + h-c--2a-R^ + dP), 

JA + 4M2 - 2nM^Nf 
c = ^ ^, (4.23) 

where we have to impose the additional constraint that the hypergeometric series ter- 
minates: 

1 + 6 - c = -AT, (4.24) 

with A^ a positive integer. To first order in the number of background D7-branes, the 
mass spectrum of the D7' probe scalar fluctuations is 

M^ = 4(A^+ |m| + |n| +2)(A^+ |m| + |n| + 1) 

+A^/(2(A^+ \m\ + \n\ + 2)(A^+ \m\ + \n\ + 1) 
(3 + 2(|m| + \n\) + 2N){n^ + \n\ + \mn\ - 1) 



7r(|m| + \n\ + 1) 



(4.25) 



Notice that by setting Nf to zero we recover the spectrum of the scalar fluctuations of 
a D7' probe in the AdSs x S^ derived in [9]. Holographically, the scalar fluctuations 
in [9] correspond to the spinless mesons of an A/" = 2 SYM gauge theory, whereas in 
our case, the mesons belong to a theory with reduced supersymmetry (A/" = 1) and 
flavor. 
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5 Review of the conifold as a Kahler manifold 

After exploring the basic D3/D7 and D3/D7/D7' systems, we now turn to the conifold. 
Here, we review the construction of the Ricci-flat metric on the conifold. We anticipate 
that the level of detail included here will provide a better understanding of the choices 
that we make in the next section when considering the inclusion of D7-branes on 
the conifold. We follow the presentation of [26] but attempt to make some of the 
computations more explicit. 

We begin by noting that the conifold is defined as the following quadric in C^: 

4 

Y^iwr = 0, (5.1) 

a=l 

where w"" are complex coordinates. It is convenient to define a radial coordinate as 

a=l 

The symmetries and general form of the conifold can be made manifest by writing the 
defining equation as: 

det W = 0, i.e. 2:1^2 — -23-24 = 0, (5.3) 

where 

2 \wi + iw2 -ws + iw^J yzi Z2J 
In this case, the radial coordinate is given by 

r2 = tr(>V^>V). (5.5) 

We now clarify the Kahler structure, with an eye toward generalizations that include 
non Ricci-fiat metrics. If a metric on the conifold is Kahler it can be written as 

gfip = d^dyT, (5.6) 

where JF is the Kahler potential. A Kahler potential invariant under SU(2) x SU(2) 
will be a function of r^ only. For such Kahler potentials, the metric can be written as 

g,, = T' d,d,r^ + roy d,r\ (5.7) 
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where prime means differentiation with respect to r^. In terms of W the metric takes 
the form 

ds"^ = J^'tr (dW^dW) + J^"\trWUW\'^. (5.8) 

The Ricci tensor for a Kahler manifold is given by 

R^,p = -df,dp\ogg, (5.9) 

where g = det (^^p. The Ricci flatness condition is easily obtained from computing the 
determinant of the metric (5.7). In particular 

g = detg,, = -\ {r\rf + r\rfT") . (5.10) 

As not all coordinates are independent, one of the coordinates, say W4, may be elimi- 
nated in favor of the other three complex coordinates, i.e. w\ = —w\ — w\ — w\. 

The question of finding the Ricci flatness condition can now be rephrased as finding 
the Kahler potential such that (5.10) is a product of a holomorphic and antiholomorphic 
function {e.g. \i^A\^). The simplest way to achieve this is by introducing 7 = r^jF'. The 
determinant can now be written as 

^ = ^^nT^(7T (5-11) 

It now follows that choosing 7 oc r^^^ will remove all r dependence. This yields the 
Kahler potential 

jr = (^2^|2/3 ^ (^^^^^ ^ ^^^^ ^ ^^^^ ^ WiWif'^. (5.12) 

After introducing p = y^r'^^^, we obtain the standard Ricci-flat metric on the 
conifold: 



A ' 1 ' A 

ds"^ = dp^ + p^ i- '^ide'f + sin 6',^rf0- ) + q (# + XI ^°^ ^i^cpif 

\ 1=1 1=1 / 



(5.13) 



The angular variables arise from the natural SU(2) x SU(2) symmetries of the conifold. 
Furthermore, they allow us to write the general solution of the defining equation (5.3) 
as 

W = rL,Z,Ll (5.14) 
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where the Lj's are SU{2) matrices and Zq is a particular solution to the defining 
relation; for example 

_ /cosi^je5(^^+'^^) -sini^je-5(V'.-'^^)\ ^ _ /^O A 

^' ~ \^sini^,et(^^-^^) cos|^,e-i('^:'+'^i)J ' ° ~ 1^0 j ' 

Note that, as a result of the coset structure of the conifold, only the combination 
"^ = "^1 + "^2 appears in the metric (5.13). 

6 D3/D7 on the conifold 

The idea of considering a stack of D3-branes at the apex of the conifold was originally 
considered in [27] as a way to obtain A/" = 1 superconformal gauge theories. This model 
is very attractive from a variety of viewpoints, and has been subsequently developed 
into one of the trademark examples of the gauge/gravity correspondence [3] . 

In this section we consider adding D7-branes to the Klebanov-Witten background. 
We explicitly develop two ways of adding D7-branes which are characterized by the 
method of embedding of the D7-branes. 

Our approach to including backreacted D7-branes in the geometry of the conifold 
is rooted in preserving the Kahler structure. The supersymmetry and integrability 
conditions give, as before, the same condition presented in (2.11). 

6.1 The t(74 = embedding 

Given the defining equation of the conifold, (5.1), a natural embedding is to place 
the D7-branes at one of the w's equal to a constant. Throughout this section we 
consider the embedding where the D7-branes are defined by w^ = 0. Note, however, 
that we can use the SO (4) symmetry of the conifold to substitute W4 for any other 
Wi- The defining equation of the conifold and the definition of its radial coordinate 
preserve the SO (4) symmetry. This symmetry prevents us from finding a simple Kahler 
potential corresponding to the inclusion of the D7-branes. In particular, it can be 
shown that there are no modifications of the Kahler potential (5.12) of the form JF = 
{r'^f{wi,W4) + h{w4,W4)y^^ for arbitrary functions / and h. 

Being confronted with this fact rules out the most natural modification to JF. Thus, 
taking into account the symmetries of the conifold, we are led to consider solutions in 
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two variables of the form T = J-'{s, t) where 

s = wiWi + W2W2 + wsivs, t= Jwl + wl + wlJw'l + wl + vol. (6-1) 

We now seek a solution to Eq. (2.11) with Imr = 1 — n/ logt, where n/ = A^//(47r): 

det{d„AnJ') = nnimr = ('^^ ^-njlogt ^ ^^ 2) 

The denominator and constant factor on the right hand side are to ensure that the 
Ricci-flat solution JF = (r^)^'^ works for n/ = 0. With the Kahler potential taken to 
be a function of both s and t, this equation becomes 

\{s'-t^)T,{T,s{Ttt^t-^Tt)-Tl) = ("0 (1-n/logt). (6.3) 
We look for a linearized solution of the form 



J^(s,t) = (s + t) 



2/3 



l+n^^(l)+nJ^(2)+^3^(3) + 



(6.4) 



Note that the zeroth order solution naturally corresponds to the conifold in the absence 
of D7-branes. 

We now observe that the linearization of (6.3) splits the perturbation problem into 
an infinite set of equations that can be schematically written as 

L^« = ^logt, 

LJ^^(") = 5„[j^W,...,J^("-i)], (6.5) 

where Si are functions that can be explicitly determined. The most important property 
is that at any order i, the function Si depends on already determined functions of lower 
order in the perturbation expansion. In addition, the same second order differential 
operator 

L = (3s^ + 2st - 2f)d^^ + (2st + t^)d^ + (2st + 4:t^)dsdt 
+ (9s + 6t)ds + {2s + 5t)dt + 4 

, „/o / \, ,,/., /3s2 + 2st-2t2 st + 2t^\ fds\ , , 

. (» + *)-'= (a. a,)(..«p(^ ^^^^^, 2si^e)\sy' '^■^' 
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appears at every level in the perturbation expansion. Given that the operator L is 
linear, and given the general form of the sources, we conclude that the above system 
is consistent. Moreover, one can, in principle, find a solution for any given set of 
boundary conditions. 

In general, the functions JF^") depend on both s and t. However, it turns out that 
for JF(^) we can find a solution of the form 

^« = i(logt-l). (6.7) 

At second order in perturbation theory the situation is radically different, as there is 
no solution depending only on t. The equation at second order is 

i^'^'^(|)*[i^-log^*]. (6,8) 

To solve the second (and higher) order equation in general, we would like to find a 
Green's function for the linear operator L. To do so, it is convenient to perform the 
following change of variables: 

s = pcos^9, t = psm^9. (6-9) 

In this case the operator L in (6.6) becomes 

L = Sp'^dl + llpdp + 4 + (1 - i sec^ 6)81 + |(cot 6 - tan ^(3 + sec^ e))de. (6.10) 

A general solution to the homogeneous equation can be constructed using separation 
of variables. Namely we take 

:F'~'\p,e) = R{p)x{e). (6.11) 

The 9 equation, Lx{6) = Ax(^), is solved by Legendre polynomials 

X = Pn{-S + 4cos'9), n= '^^~\ (6.12) 

Note that, because of the way s and t are defined in (6.1) the angular variable is 
restricted to | < cos^ 9 < 1. Thus the argument of the Legendre polynomial lies in the 
range [—1, 1]. The remaining R{p) equation is then homogeneous in powers of p; hence 
it may be solved by taking R r^ p^. This gives the general solution to the homogeneous 
equation of the form: 



oo 



^(2) = y c±p-|(2±Vi+3«(«+i))p^(_3 + 4 cos^ 9), (6.13) 



n=0 
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where the c^ are coefficients that must be determined by boundary conditions. Of 
course, to obtain the actual solution to the inhomogeneous equation (6.8), we would 
have to proceed with an actual construction of the Green's function. Nevertheless, even 
without doing so, we can discern some general features of the solution. For example, the 
role of n in (6.13) can be understood as follows. Given that p = s + t is the distance 
from the D3-brane, it is natural to conclude that higher terms in n correspond to 
deforming the gauge theory by operators with higher conformal dimension. 

From the Kahler potential we can easily derive the metric to ffist order in the num- 
ber of D7-branes, Nf. "Turning on" the D7-branes alters the g^^ conifold geometry 
in two ways: 

9mn = gin (l + ^(logt - 1)) - ^r'''drnw'dnw\ (6.14) 

First, there is the rescaling of the background metric with a factor (l + A^/(logt — 1)/3). 
Second, we notice the presence of some off-diagonal metric components due to the 
last term in (6.14). More explicitly, to ffist order in Nf, the correction to the six- 
dimensional conifold geometry reads 



N. 
^ 



/^4/3 



-(1 — cos^i COS ^2 ~ 2 sin^i sin ^2 cos(0i + (p2))dr'^ 



-\ — r(sin 62 cos 9i — cos 61 sin 62 cos(0i + (j)2))dr dOi 

+ -r(sin 9i cos 62 — cos 62 sin 9i cos(0i + (j)2))dr d92 

+ -r sin 61 sin 62 sin(0i + (f)2)dr d{(f)i + ^2) 
rV 

+ — ( (cos 61 - COS 92)dip d{(j)i + 02) 

+ sin 62 sin(0i + (t)2)dip dOi — sin 61 sin(0i + (p2)d'ip ^6*2 

+ sin 9i cos 02 sin(0i + 02)c?^2 d{(f)i + ^2) 

-|- sin 62 cos 61 sm{(f)i + 02)c^^i d{(f)i + ^2) 

+ (1 — cos 9i cos 02 + sin 9i sin 62 cos(0i + (j)2))d{<t>i + ^2)^ 



(6.15) 



where w^ is given by (7.13). 
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6.2 A P^-inspired holomorphic embedding 

Here we consider a different embedding which is inspired by the small resolution of the 
conifold [26,28]. In this embedding the D7-branes are perpendicular to a P^ subman- 
ifold inside the conifold. The resolution of the conifold can naturally be described in 
terms of the {zi, Z2, z^, z^) of (5.4). In particular, resolving the conifold means replacing 
the equation Z1Z2 — Z3Z4 = by the pair of equations 

= 0, (6.16) 

where the product A1A2 7^ 0. Note that (Ai, A2) G CP^ (any pair obtained from a given 
one by multiplication by a nonzero complex number is also a solution). Thus (Ai, A2) 
is uniquely characterized by the ratio A = A2/A1 in the region where Ai 7^ 0. Working 
on this patch, a solution to (6.16) takes the form^ 

W=(~''^ '']. (6.17) 

\-Z2\ Z2J 

Thus (z2,-23,A) are the three complex coordinates characterizing the resolved conifold 
in the patch where Ai 7^ 0. 

As usual, the conifold metric is Qmn = dmdnJ^, where JF is the Kahler potential. In 
this case, the radial coordinate defined by (5.5) takes the form 

r^ = tr(W^W) = (1 + |An(|^2r + ksH- (6.18) 

The Kahler potential for the conifold is simply^ 



jr = ^^2^2/3 



(l + |A|^)(|^2|^ + k3|^) 



2/3 



(6.19) 



In this parametrization it is natural to consider a modification of the Kahler potential 
that depends on the two variables: 

s = \Z2\^ + \z^\\ t = \X\\ (6.20) 



'*In the region where Ai is aUowed to be zero we have instead A2 7^ 0. In this case the general 

solution can be written as W = , where ^ — A1/A2. 

^Note that in this specific patch the matrix (6.17) satisfies the defining equation for the conifold. 
One can think of this matrix as a "solution" with zi = — AZ3 and 24 = — Az2- 
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Doing so, we find that the determinant of the metric takes the form 

det gmn = i^t^s + ^t^s^ssS - ^s^lt^s + TttJ^lt + J^s^ss^ttts) . (6.21) 

As explained previously, our task is to find a solution J-'{s,t) of equation (2.11) with 
the condition that when the number of D7-branes goes to zero we recover the Kahler 
potential of the conifold. 

We now observe that the above equation is homogeneous in powers of s (as well as 
in powers oft). Hence, although the equation is nonlinear, in principle all s dependence 
could be removed by a judicious choice of JF ~ s" for some power n. 

This suggests a separation of variables in s and t of the form JF = s"'f{t). Eq. (6.21) 
becomes: 

detgmn = nh^'^-'^ f{ff"t + ff - {f')H). (6.22) 

Taking n = 2/3 in the above gives an s independent determinant, and leaves a non- 
linear ODE to solve: 

detgmn=(^l^ /(//"t + //' - (/O'^) = Q) (1 - n^log(t)). (6.23) 

We have fixed f^O = (2/3) by demanding that / reduces to /(t) = (1 + t)^/^ for Uf = 
{i.e. in the absence of D7-branes). Again, we have assumed that ^2 does not change 
when Uf is nonzero. Although the source does not respect homogeneity in t, it is still 
useful to perform the change of variables 

t = exp{y), fit) = expiy/3)l{y), (6.24) 

In this case, (6.23) implies that 

l%dylog{l) = {2/3){l-nfy). (6.25) 

The standard conifold case can be recovered by taking / = (2 cosh(|))^/^. Interestingly, 
the combination / = (2Acosh(^^^))^/'^ also gives a Ricci fiat metric where the param- 
eters are related to those discussed at length in [29] ; however we take the simplest case 
with A = 1, yo = 0. 

Note that, in fact, the l^ prefactor in (6.25) may be replaced by any non-zero power 
of / through the replacement / — > /" where a is an arbitrary non-zero constant. Of 
particular interest is to take / -^ /^/^, because the conifold itself, (6.19), has this 2/3 
power built in. This gives the equation 

fdydy\ogil) = {l-nfy). (6.26) 

27 



The Ricci-flat solutions (corresponding to n/ = 0) are now simply / = 2cosh(|). 
Assuming a perturbative solution to (6.26) of the form 

7/ 

/ = 2 cosh(-) exp{nfgi{y) + n)g2{y) + ■■■), (6-27) 

we find that at each order, we must solve an equation of the form 

Lgn{y) = Sn[gi,...,gn-il L = 1 + (1 + coshi/)9j. (6.28) 

Noting that L admits the homogeneous solution 

g{y) = ci tanh(|) + c^ (l/tanh(|) - 2) , (6.29) 

we may construct the Green's function 

(2 - yo tanh f ) tanh | - (2 - |/ tanh |) tanh f ^^ 

G(l/,l/o) =^o(l/) + ^ -^ — r ^ -0{y-yo), (6.30) 

1 + cosh I/O 

where (70(1/) satisfies the homogeneous equation, and must be chosen to respect bound- 
ary conditions. 

At linear order in n/, we find Lgi = —y/2, which admits a simple solution 

91 = -f- (6.31) 

Taking this form of gi, the next order equation has source Lg2 = —y'^/A, and may be 
solved to yield 

g2 = -2 log(l + e^) + ^, ^ + - — -Li2(-e^), (6.32) 

^ sv ; 2(1 + 6?') 1 + ey ^ ^ ^ ^ 

up to the possible addition of a solution to the homogeneous equation. 

We have therefore managed, in this simple case, to find a modification to the 
Kahler potential of the conifold accommodating the presence of D7-branes up to second 
order (and in principle extendible to arbitrary order) in the number of D7-branes. 
Interestingly, however, as we will explain in the next section, this embedding of the 
D7-branes seems to be in a different class from the W4 = one. 

7 Gauge theory 

Let us start by commenting on the decoupling limit. To make connection with gauge 
theories we first show that the supergravity solutions presented here allow for a region 
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where only the low energy D-brane degrees of freedom are described. This is a necessary 
requirement to claim a duality between our supergravity backgrounds and the gauge 
theories appearing in the worldvolumes of branes. 

The Yang-Mills coupling of the gauge theory living on a Dp-brane is 

9lp = 9sir'- (7.1) 

We want to take the string scale to zero Is ^ keeping the D3 SYM coupling fixed 
9d3 ~ 9s- This automatically means that g^j~,-j -^ which implies that the D7-brane 
dynamics decouple, leaving us with a flavor symmetry. 

There are other quantities that we will keep fixed in the limit /^ — > such as the 
mass of the 7-7, 3-7 and 3-3 strings. This translates into a rescaling of the radial 
coordinate perpendicular to the D3-brane as u = r/l'^. This decoupling limit provides 
us with a consistency check on the solutions we find. In particular, we will require that 
when turning off the number of D7-branes the metric becomes either AdSr, x S^ for 
the D7 in flat space or AdSr, x T^'^ for D7's on the conifold. 

The gauge theory of a stack of D7's on a D3 background has been discussed in 
[21]. A complete account of the degrees of freedom of this system (and in general of 
Dp/D{p+4)) is given in the textbook [30]. One interesting fact about the D3/D7 system 
is that the condition for the dilaton to be small, that is for the string loop corrections 
to be small requires r <^ ^'^^KasNf) ^ The condition to suppress the a' corrections, on 
the other hand, implies in the case of D3/D7 the radial distance from the D7 must 
satisfy r ^ ^-^nNc/Nf ^ These two conditions have to be taken into consideration when 
interpreting the range of validity of the supergravity calculation. Generically, for large 
gsNf we must consider only Nc ^ Nf. However, for very small gsNf, there is a range 
in which, in principle we could relax the relation between Nc and Nf. 

7.1 Matching Pathologies: Naked singularities and Landau Poles 

If one restricts to the form of the dilaton given in 2.16 and rewrites the metric using 

( Nf 



(isg = dzidzi + dz2dz2 + I 1 — -. — log -232:3 dz^dzs, (7.2) 

it becomes clear that the metric has a naked singularity when 



1-^ln^ , (7.3) 
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vanishes, that is, at: 

PL = Poe2-/^=^/. (7.4) 

A simple way to reahze the relation to the Landau pole is through the identification 

^ e-*. (7.5) 



9ym 



Thus, the metric singularity coincides with the point at which qym — *■ cxd. It is, 
therefore, natural to consider only the region with ps < Pl- 

Gauge theory dual of D3/D7/D7' 

The gauge theory living on this brane configuration has been described in the literature, 
for example in [24,4]. The analysis is closely related to a configuration of intersecting 
M5's discussed in [31]. The theory includes the A/" = 4 supermultiplet plus Nf hyper- 
multiplets Q and Ni hypermultiplets T transforming in the fundamental representation 
of SU(A^). A simple way to describe this theory is via the superpotential 

W = Wm=a + \qQZ2Q + AtTZsT, (7.6) 

where Z2 and Z^ are two of the three chiral superfields of the A/" = 4 multiplet and we 
are explicitly using A/" = 1 notation. 

The results of our supergravity analysis point to a generic appearance of nondiag- 
onal terms in the metric. The explicit form of the metric (3.16) shows a complicated 
pattern of mixing the Z2 and the z-^ coordinates. This seems to suggest that, at strong 
coupling, the two fiavor sectors mix in a complicated way due to interactions. 

Gauge theory dual of D3/D7 on the conifold 

D3-branes on the conifold [27] realize an A/" = 1 superconformal field theory (eight 
supercharges) with gauge group SU(N) x SU(N). The matter content of this theory 
consists of chiral superfields Ai and B-i with i = 1,2 transforming as (N, N) and (N, N) 
respectively. The fields A and B are doublets under one of the SU(2)'s in the global 
symmetry SU(2) x SU(2). Under this symmetry a natural superpotential to add is [27] 

W,^,foid = ^e^'e'^^TrABkA.Bi. (7.7) 

As the example of A/" = 4 shows, the structure of the deformation of the superpo- 
tential in the gauge theory is dictated by the geometry in the supergravity background. 
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In particular, note that the superpotential for the conifold is intimately related to the 
defining equation for the conifold (5.3) under the identification 

zi = AiBi, Z2 = A2B2, Z3 = A1B2, Zi = A2B1. (7.8) 

To write down the superpotential corresponding to the embedding defined by W4 = 0, 
recall that 1^4 ~ 2:1 + Z2 and therefore the natural deformation of the above superpo- 
tential is of the form 

W = Wconifoid + XqQ{AiBi + A2B2)Q + Xgq{BiAi + ^2^2)9. (7.9) 

In [7] a deformation of this type was motivated using the theory of D3-branes on C^/Z2. 

Comments on the /3-functions 

The exact [3 function for the canonical gauge coupling of A/" = 1 with gauge group 
SU{N) and F flavors is [32] 

lovr^ 1 _ NJ-^ 

where 7 is the anomalous dimension of the matter fields. One of the most attractive 
features of the solutions describing branes on the apex of the conifold is the precise 
geometrical description of the NSVZ /9-function [27,33]. 

One very important feature of (7.10) is that it serves as a guide in determining the 
boundary conditions of the supergravity solution describing fiavor. Namely, as describe 
by equation (2.14) the dilaton can be any holomorphic function of the coordinates. In 
particular the sign of the coefficient of the term containing the logarithm is arbitrary 
in the supergravity solution. In this case equation (7.10) dictates the sign. 

Another important remark about (7.10) is that it also allows to understand the 
region of the radial coordinates where the theory is describing the fiavor sector. Note 
that in general the dilaton goes from a constant at infinity to a logarithmic term 
near the D7-branes. As we will see in what follows for the case of the conifold, only 
the logarithmic region (decompactification limit) contributes to the /3-function in the 
expected way. 

For the Klebanov-Witten theory the assignment of anomalous dimension of 7 = 
— 1/2 to the operators TiAiBj guarantees that the beta function is zero, which corre- 
sponds to a superconformal fixed point [27] 
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^^3iV-2iV(l-7). (7.1i: 



o?log(A//x) gj 
The two gauge couplings are related to the supergravity background as [27,3]: 

47r^ Att'^ it 



-e-\ 



gl gl gs 



Att^ A-n^ 1 



gi gi Svra' 




n mod 271. (7.12) 



The case with probe D7-branes was considered in [7]. In the absence of B2 field we 
obtain that the two couplings are the same. To read off their value we need the explicit 
expression for w^ which can be obtained from (5.3), (5.14) and (5.15) 

^4 = ILef (^+<^i+<^^) [cos ^ sin ^ - sin ^ cos ^g-^^^^+^^n . (7.13) 

v^ L 2 2 2 2 ^ ^ ' 

Assume, that the renormalization energy scale is proportional to the radius A ~ p (the 
perpendicular distance from the D3-branes) we have the following expression for the 
running of the coupling 

-^ = 4^/(1 -7g), (7.14) 



d^g^M 2 

The value of 7q can be found, similarly to a computation of [7], by power counting. 
Given that the dimension of the superpotential (7.9) is three, and that the dimension 
of the fields A's and 5's is 3/4 as argued in [27], we conclude that 7q = 1/4. 

For the case of the P^ parametrization of the conifold we note that r ~ log A and 

A = tan— e-*'^^ (7.15) 

This combination is independent of the radial coordinate p, therefore, keeping the 
natural identification of p with the RG energy scale A we conclude that this deformation 
of the gauge theory does not affect the /9-function and therefore does not truly amount 
to adding flavor to the KW gauge theory. 

8 Conclusions 

We conclude with a list of open problems which we believe would be interesting to 
address. The first, and perhaps most haunting problem in the gauge/gravity corre- 
spondence is how to achieve asymptotic freedom, or in other words how to follow the 
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theory from strong coupling where supergravity is a good approximation into the weak 
couphng regime. This problem is well beyond the scope of the present paper. But it 
must be kept in mind in order to understand the limitations of the present work. Some 
direct open problems resulting from our investigation are: 

• Application of our techniques to other Kahler spaces. Most prominently, the 
resolved and deformed conifolds. 

• The inclusion of 3-form fluxes. In the case of the conifold we refer to G3 which 
corresponds to including fractional branes which change the relative rank of the two 
gauge groups. 

• An interesting open problem would be to understand the backreaction for theories 
that do not necessarily admit a formulation exploiting Kahler geometry. Such seems 
to be the case for the Maldacena-Niiiiez solution and the D4/D6 bound state. 

• Ideally, it would be desirable to obtain the spectrum of the different modes for 
the class of solutions we have presented. We should point out that some analysis of 
the spectrum (for example of mesons) has been presented in the probe approximation 
[11,9,10,8]. But an analysis of the fully backreacted case will certainly yield a more 
robust picture. In this paper we took a step further by showing how the backreacted 
D3/D7 solution could be used to derive the scalar (meson) spectrum of a A/" = 1 gauge 
theory realized on the worldvolume of a D7' probe. 

To conclude, we believe we have explicitly demonstrated a viable way to construct 
fully backreacted supergravity backgrounds dual to four-dimensional M = 1 SYM with 
flavor in the fundamental representation. We hope to return to many of the fascinating 
issues in the gauge/gravity correspondence with flavor in future work. 
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A Comments on the Monge- Ampere equation 

In this appendix we collect some known results about the Monge-Ampere equation 
that are relevant for the analysis in the main text. The literature on this subject is 
vast [34-36] and we refer the reader to those references for further details. 

We start with some motivation of why Monge-Ampere equations are relevant for our 
analysis. We assume (four-dimensional) Poincare symmetry and (for simplicity) turn 
off 3-form flux. As in section 2, one arrives at equation (2.11) for type IIB supergravity 

det(a^apJ^) = finim(r). (A.l) 

To place this in the proper context of the Monge-Ampere equation, we must change 
the differentiation appearing above into differentiation with respect to real variables. 
First, we consider the "decompactification" limit of seven branes [19], where r ~ 
iC + ialog{z). We will then place our D3's close enough to the D7's so that this 
approximation is valid. Thus Im(r) ~ log{zz) where we have chosen our origin so 
that 6 = 0. We then make the Ansatz that the functional dependence of the Kahler 
potential is J-' (\og{ziZi), . . . ,\og{znZn)), where the log's have been inserted for later 
convenience. This requires that Q has the form Q = Ilz^^ for any set of m's so that 
both left and right hand side are functions of the new variables Un = log(z„z„). Under 
this change Eq. (A.l) becomes: 

det((9j(9jjF) 
nexp(i/i) 
det{didjJ^) = exp(S(mi + l)|/i)Im(r). (A.2) 

Note that the change of variables has introduced only a product of the z^s to powers. 
This can be absorbed into the definition of Q, and so has no affect on the curvature. 
Equation (A.2) is the standard Monge-Ampere equation in real variables. We now turn 
to some of the main results from the study of Monge-Ampere equations. 

The Dirichlet problem for Monge-Ampere equations is formulated as follows. Let 
r2 be a bounded strictly convex domain in R" {n > 2) defined by a C°° strictly convex 
function h on Cl satisfying h\gn = 0. Given u{x), a C°° function on dQ which is the 
restriction to dQ of a C°° function 7 on Cl, we consider the equation: 

logdet(9i(9j0) = /(x,0), (pldn = u, (A. 3) 

where f{x,t) G C°°(r2,]R). The main result that we will use is an existence theorem 
which can be stated as follows [34-36]. The Dirichlet problem (A. 3) has a unique 
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solution belonging to C^{Q), when n = 2, if there exists a strictly convex upper 
solution 7o G 0^(0,) satisfying 

det((9j7o) < exp/(x,7o), 7o|9n = u, (A.4) 

and if fi{x, t) > for all x G ^2 and t < supg^ u. 
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